Partial-wave contributions to pairing in nuclei 



Simone Baroni, 1,2 'E] Augusto 0. Macchiavelli, 3, [l| and Achim Schwenk 2 ' 4 ' 5, [l 

1 Institute for Nuclear Theory, University of Washington, Box 351550, Seattle, WA 98195, USA 

2 TRIUMF, 4004 Wesbrook Mall, Vancouver BC, V6T 2A3, Canada 
s Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum fur Schwerionenforschung GmbH, 64291 Darmstadt, Germany 
J Institut fur Kernphysik, Technische Universitdt Darmstadt, 64289 Darmstadt, Germany 

We present a detailed study of partial-wave contributions of nuclear forces to pairing in nuclei. 
For T = 1, J = pairing, partial waves beyond the standard 1 So channel play an interesting role 
for the pair formation in nuclei. The additional contributions are dominated by the repulsive 3 Pi 
partial wave. Their effects, and generally spin-triplet nuclear forces between paired nucleons, are 
influenced by the interplay of spin-orbit partners. We explore the impact of including partial waves 
beyond the 1 So channel on neutron-neutron pairing gaps in semi-magic isotopic chains. In addition, 
we show that nuclear forces favor T = 1, J = over T = 0, J = 1 pairing, except in low-j orbitals. 
This is in contrast to the free-space motivation that suggests the formation of deuteron-like T — 
pairs in N = Z nuclei. The suppression of T — pairing is because the 3 Si strength is distributed 
on spin-orbit partners and because of the effects of the repulsive Pi channel and of D waves. 
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I. INTRODUCTION 



and total angular momentum J. 



Phcnomcnological energy density functionals are im- 
pressively successful in the description of medium-mass 
and heavy nuclei but lack a microscopic connection 
to nuclear forces and seem to have reached the limits 
of improvement in the present functional form 0, 
This has lead to exciting efforts, largely driven by ef- 
fective field theory (EFT) ideas, to develop a universal 
energy density functional based on microscopic interac- 
tions 

H IS U 0, H- These developments rely on the 
Hartree-Fock approximation as a good starting point for 
nuclei or on a perturbative expansion about nuclear mat- 
ter. Both have become possible by evolving nuclear forces 
to low-momentum interactions using the renormalization 

group (rg) 0, m m m, m m . 

In this paper, we follow Refs. [lHGil and use existing 
energy density functionals in the particle-hole channel, 
to build a reasonable self-consistent single-particle basis, 
combined with low- momentum interactions Vi OVJ t @, H3] 
in the pairing channel. The calculational details are given 
in Sect. [TTJ This was shown to provide a good starting 
description of neutron-neutron and proton-proton pair- 
ing properties in semi-magic nuclei lj| [H, [l?], EH- In 
these and other [ijj [2(| calculations based on nuclear 
forces, the pairing interaction was restricted to the 1 So 
partial wave. This is motivated by pairing in infinite mat- 
ter (see, for example Ref. [U), but in nuclei the paired 
nucleons are not in back-to-back-momentum configura- 
tions. Therefore, other partial waves can contribute to 
the pairing interaction for two particles with isospin T 
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In Sect. IIII1 we discuss partial- wave contributions to 
pairing interactions in nuclei. This is followed in Sect. II Vl 
by a detailed analysis of T = 1, J = pairing matrix el- 
ements in the sd and pf shells, where we focus on the 
effects of partial waves beyond the standard 1 So chan- 
nel. Spin-triplet partial-wave contributions to pairing in 
time-reversed states are found to depend on the spin- 
orbit configurations involved. In Sect.fVl we explore the 
impact on T = 1 pairing gaps in semi-magic isotopic 
chains. Our analysis of pairing matrix elements provides 
a simple explanation of the results observed in nuclei. 

We extend the study of partial-wave contributions to 
T = pairing in Sect. IVI1 and compare the pairing 
strengths to the T — 1 , J = channel at the level 
of low-momentum interactions. Since the deuteron is 
bound, the free-space motivation suggests the formation 
of T = 0, J = 1 pairs in all N = Z nuclei [II]. In 
contrast, there is no direct observation of T = pair- 
ing to date. While we do not perform a calculation for 
nuclei, our analysis of T = 0, J = 1 pairing matrix ele- 
ments based on nuclear forces can provide a microscopic 
explanation for the suppression of T = pairing. For 
the 3 Si channel, the nuclear force strength is distributed 
on the spin-orbit partners, leading to a geometrical sup- 
pression in higher-j orbitals. Moreover, the additional 
partial-wave contributions are dominated by the repul- 
sive 1 Pi channel with some smaller contributions due to 
D waves. As a result, in the absence of many-body ef- 
fects on pairing, we find that low-momentum interactions 
favor T = 1, J = over T = Q, J = 1 pairing, except in 
light N — Z nuclei where pairing is dominated by low-j 
orbitals. The standard motivation for T = pairing in 
nuclei is therefore at best incomplete. We conclude and 
give an outlook in Sect. I VIII 
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II. CALCULATIONAL DETAILS 

The minimization of the energy density functional in 
the presence of pairing leads to solving the self-consistent 
Hartree-Fock-Bogoliubov (HFB) equations [l|, [23[ that 
determine the quasiparticle (q) basis, 



h- [i A 
-A* -(h-fi) 



U" 



= E n 



U" 
V" 



(1) 



Here h denotes the single-particle Hamiltonian, \x the 
Fermi level, E q is the quasiparticle energy and U q ,V q 
the corresponding coefficients of the Bogoliubov trans- 
formation from single-particle to quas iparticle states. We 
use the Skyrme functional SLy4 |24| and first solve the 
Hartree-Fock (HF) equations, h(f> a — e^ F cf)a, on a spheri- 
cal mesh with 0.1 fm spacing and 16.0 fm box radius. Our 
results are stable with respect to increasing the radius 
and decreasing the mesh spacing. This defines the single- 
particle basis |a), using the shorthand label a = n a l a j a 
with radial quantum number n a , orbital angular momen- 
tum l a and total angular momentum j a = l a ± 1/2. 

Using the HF single-particle Hamiltonian, we then 
solve the HFB equations, where the state-dependent gap 
matrix A for T — 1 , J = pairing is given by 



'2j e + l y/(l + M(l + M 
2 ja + 1 2 



cd N q 

x (ab\(l- P 12 )Vi owk \cd}j =0 ,T=i,T z =-i- (2) 

This defines the gap equation and we focus on neutron- 
neutron (T z — —1) pairing properties. The matrix 
elements of the pairing interaction in the second line 
of Eq. ([5]) are antisymmetrized using the exchange op- 
erator P\2 and normalized, and d ab is shorthand for 
&n a n b Si a i b Sj a j b . In general, after the HFB equations are 
solved self-consistently, one has to insert the resulting 
densities back into the single-particle Hamiltonian, which 
leads to a new set of single-particle states, in turn to a 
new HFB solution, and these iterations have to be re- 
peated to obtain the fully self-consistent HFB solution. 
However, as shown in Sect.fV] the feedback on the single- 
particle Hamiltonian can be neglected for pairing prop- 
erties to a good approximation. 

For the neutron-neutron pairing interaction, we start 
from the chiral N 3 LO two-nucleon (NN) potential (A = 
500 MeV) of Ref. [I|| and use the RG to evolve this 
NN potential to low-momentum interactions Vj ow k with 
a smooth n eX p = 4 regulator with A = 1.8 — 2.8 fm~ 1 [l(| ■ 
This evolution renders the many-body calculation more 
controlled [H, [H, [H| an d provides a good starting 
point for connecting energy density functionals to nuclear 
forces [1, [H| . Based on the universality of V\ ow k [3 EH , 
we do not expect large differences starting from different 
N 3 LO potentials. As discussed in the following section, 
we calculate the jj-coupled pairing matrix elements en- 
tering the gap equation, Eq. by expanding the HF 
single-particle states on the harmonic oscillator (HO) 



basis. For pairing properties around the Fermi level, 
in particular for the lowest-quasiparticle-energy canon- 
ical state and average gaps presented in Sect. [V] we 
have found that the HFB single-particle space can be re- 
stricted to states below 60 MeV (for hard potentials with 
large cutoffs, states up to ~ 1 GeV are needed for conver- 
gence [19ll20lp. Finally, an important direction for future 
work is to use energy density functionals in the particle- 
hole channel that are based on the same nuclear forces 
and to include many-body contributions to the pairing 
interaction. First results in this direction have been pre- 
sented in Ref. [H]. 

In HFB theory, the interacting particles pair in or- 
bitals related by time-reversal symmetry [27|], where the 
radial quantum numbers can differ due to the Bogoli- 
ubov transformation [23| . Semiclassically, this is realized 
by two particles moving in the same orbit with opposite 
time order, so that (l a ,j a ) = (h,jb) and (l c ,j c ) = (k,jd) 
and hence the pairs have positive parity. For complete- 
ness, the gap equation is in the so-called phase con- 
vention II [271 ] . for which the operation of the time- 
reversal operator T on a single-particle state is given 
by T\l a j a m a ) = (-l)-?«~ m » \l a j a - m a ), with magnetic 
quantum number m a . 



III. PARTIAL- WAVE CONTRIBUTIONS TO 
PAIRING INTERACTIONS IN NUCLEI 

The jj-coupled pairing matrix elements in the HF ba- 
sis are obtained from the interaction matrix elements of 
Mowfc in two-particle spherical HO states, which are cal- 
culated from the different partial-wave contributions us- 
ing the standard recoupling formula, 



(a b | (1 - P 12 ) Ui ow fc | cd)j,T,T„ 
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(3) 



Here we use x = y2x+l and standard notation for 
6j and 9j symbols [28j. The bracket (nlN L\n a l a nblbX) 
denotes Talmi-Moshinsky brackets [29[, with relative 
and center-of-mass radial and orbital quantum numbers 
n, I, N, L, and total orbital angular momentum of the pair 
A. The two-body spin is given by S and J Te \ = I + S de- 
notes the relative total angular momentum. We use the 
standard notation 2S+1 l j rol to denote the different partial 
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waves, which we list for completeness: 
forT=l: % 3 Po,i,2 X D 2 3 F 2 , 3 , 4 l G 4 



forT = 0: 3 Si 1 P 1 3 B 1 



2,3 



3 G 



3,4,5 



Additional restrictions come into play in the pairing 
problem when the total angular momentum of the paired 
nucleons is J — 0. Because L + J rc \ = J, one therefore 
has L = J rc i for J — 0. Moreover, because the parity 
of the pair is positive, it follows that the relative and 
center-of-mass orbital angular momentum I and L must 
have the same parity. As a result, the relative orbital and 
total angular momentum I and J le \ also have the same 
parity, and therefore only uncoupled channels 2S+1 Zj rol= / 
contribute to T — 1, J — pairing matrix elements. On 
the other hand, T = 0, J — 1 matrix elements are not af- 
fected by these additional constraints and coupled chan- 
nels also contribute to the pair formation in nuclei. 

In summary, this shows that the 1 So and 3 Si — 3 Di 
partial waves are not the only channels contributing to 
T = 1 and T — pairing in nuclei. In addition, the 
following partial waves are part of the pairing interaction: 



T=1,J = 0: % 3 Pi l D 2 
T = 0,J=1: 3 Si 1 P 1 3 D lj2 , 3 



3,4,5 



Interestingly, both 3 Pi and Pi partial waves are repul- 
sive at the relevant energies [301 - W e therefore expect a 
reduction of the pairing gap compared to studies based 
on the standard S-wavc interactions. 



IV. T = 1, J = PAIRING MATRIX ELEMENTS 

The sd and pf shells are very useful for understand- 
ing the results observed in nuclei in Sect. |V] These two 
major shells are small enough to allow a simple study 
and at the same time they are complex enough to draw 
general conclusions. In Figs. [T] and [51 we show the jj- 
coupled T = 1 , J = pairing matrix elements in the sd 
and pf shells based on the smooth V\ ow k evolved from 
the N 3 LO potential of Ref. 25] to A = 2.0 fm" 1 . Here, 
we give neutron-proton (T z — 0) matrix elements in or- 
der to compare to the T = 0, J = 1 pairing strengths in 
Sect. IVI1 The contributions from isospin-dependent nu- 
clear forces that distinguish between the T = 1,T Z = 
and T = 1,T Z = —1 channel (used in the calculations for 
nuclei) are small and do not change this analysis. The 
matrix elements are in the HO basis with hut — 10 MeV, 
but the general effects of partial waves are qualitatively 
similar for other values of fuv or in the HF basis. 

As expected based on NN phase shifts, the pairing ma- 
trix elements are attractive for the channel. In ad- 
dition, we find that the attractive contribution in- 
creases approximately linearly with increasing j a or % 
within each (l a ,h) group. This can be understood using 
a semiclassical picture. In the classical limit, the plane 




FIG. 1: T = 1, J = pairing matrix elements in the sd 
shell. Results are shown for the RG-evolved N 3 LO potential 
of Ref. [25] with A = 2.0 fin -1 , represented in two-particle 
HO states with Yiuj = lOMeV. The label (LjJ 2 ~ (hj b ) 2 
denotes the bra— ket quantum numbers and we have grouped 
matrix elements according to (l a , h)- In the sd shell, the two- 
particle quantum numbers combined with J = restrict the 
summation over partial waves to I ^ 2. 
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FIG. 2: Same as Fig. [T] but for the pf shell, where partial 
waves with / ^ 3 can contribute to T = 1, J = pairing 
matrix elements. 



of motion of a particle is determined by the orbital an- 
gular momentum vector £. Therefore, two particles mov- 
ing in time-reversed orbits with total angular momentum 
zero will have opposite angular momentum vectors. The 
larger the value of l ai h, the closer the situation will be 
to the classical picture and the more aligned the orbital 
angular momentum vectors will be (see also Ref. [3l|, 
Sect. 2.2]). Hence, the overlap between the orbitals is 
larger with increasing orbital angular momentum. The 
same holds with increasing total angular momentum and 
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FIG. 3: Illustration of two particles in time-reversed states 
with j> —1 + 1/2 (left) scattering to the same (j>) 2 configu- 
ration (top right) or to the (j<) 2 configuration involving the 
spin-orbit partner j < = I — 1/2 (bottom right). 



the larger overlap leads to the more attractive con- 
tribution for larger j a or jb within each (l a , lb) group. 




N 



FIG. 4: Cutoff dependence of the LCS gaps in the tin isotopes 
including only the So contribution to the pairing interaction. 
The order in the legend corresponds to the order of the curves. 
The crosses denote the results of Ref. |l6j ] , based on a smooth 
Viow k evolved from the Argonne vig potential to A = 2.5 fm -1 



with n c 
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A. Spin-triplet contributions 

The 3 Pi partial wave is repulsive at the relevant en- 
ergies 30]. However, the contributions to the pairing 
matrix elements in Figs. [T] and [2] are clearly more com- 
plicated. While the 3 Pi contribution is repulsive be- 
tween the same spin-orbit partners (j>) 2 — (j>) 2 and 
(j<) 2 — (j<) 2 (as expected from phase shifts), we find 
an attractive contribution to the pairing matrix elements 
connecting {j>) 2 ~{j<) 2 and (j<) 2 ~(j>) 2 configurations, 
where the spin-orbit partners are denoted by j> = I + 1 /2 
and j< = I — 1/2. This effect can be seen prominently in 
the pf shell in Fig. [2] This is a general property of spin- 
triplet nuclear forces and applies also to T = 0, J = 1 
pairing matrix elements. 

The different configurations for scattering of two par- 
ticles in time-reversed states are illustrated in Fig. [3] 
In the semiclassical picture, the spin s of the particle 
will be in the direction or opposite to the orbital angu- 
lar momentum vector t for the j> or j < configuration, 
respectively. If the two particles scatter from the j > or- 
bital to its spin-orbit partner j < (left to bottom right in 
Fig. [3j), they will experience a spin flip that would not 
occur in the case that the particles remain in the same 
j > orbital (left to top right in Fig. [3]). This spin flip 
leads to a relative sign in the spin part of the wave func- 
tion |TT(T)) =|U) = (1-5 = 1) +\S = 0))/y/2 (for the 
left configuration in Fig. [3]) versus ||T(|)) = — |J,f) = 
(-\S = 1) + \S = 0)) /V2 (for the bottom right configura- 
tion), which results in the relative sign difference of spin- 
triplet partial wave contributions to (j>) 2 — (j<) 2 and 
{]<? ~ U>) 2 versus (j>) 2 - (j>) 2 and (j<) 2 - (j<) 2 con- 
figurations. This also holds for (l a j a ) 2 — {hj b ) 2 matrix el- 
ements with l a ^ lb, leading to an additional sign for gen- 



eral upper-lower spin-orbit configurations {j a >) 2 ~ C?6<) 2 
and (j a< ) 2 - (j fc >) 2 . 



V. T — 1 PAIRING GAPS 

As measures of T = 1 pairing in nuclei, we calcu- 
late both the lowest-quasiparticle-energy canonical state 
(LCS) gap and an average gap. The former is defined as 
the diagonal pairing gap A aa of Eq. © , where the LCS 
state a is the canonical state with the lowest quasiparticle 
energy E a given by 

E a = V(e a - m) 2 + AL . (4) 

The canonical basis is defined by diagonalizing the den- 
sity matrix, and e a is the canonical single-particle energy. 
The LCS state is close to the Fermi level using the defi- 
nition of Eq. (0| . We introduce the average gap as 



Ea 




{2j a + 1) A aa 




(2ja + 1) 



The factor (J2 q U^V^) (2j a + 1) weighs states according 
to the degeneracy of the level and averages over an energy 
window around the Fermi level, which is approximately 
given by one shell above and below the Fermi level (for 
example, « /z ± lOMeV in 120 Sn). 

We compare the theoretical LCS and average gaps cal- 
culated for even-even nuclei to experimental gaps based 
on the three-point mass formula centered on odd nuclei 




FIG. 5: Comparison of the LCS gaps in the tin isotopes ob- 
tained using a smooth Viowfc evolved from the Argonne via 
potential to A = 2.5 fm -1 with n oxp = 6 and those of Ref. [l6l ] 
using the same pairing interaction. Experimental gaps based 
on the three-point mass formula are shown for comparison. 
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using experimental binding energies B(N, Z) of Ref. [3J 

For low-momentum interactions, superfluid properties 
are dominated by states around the Fermi level and the 
contributions from the right-hand side of the gap equa- 
tion, Eq. (J2J, fall off as particle pairs scatter to higher- 
lying states [13]. Therefore, our results are independent 
of the HO basis parameters, as long the HF states around 
the Fermi level are reproduced. We have found that 10 
oscillator shells are sufficient for 6 < Two < 18 MeV. The 
dependence of the pairing gaps on the HO basis param- 
eters is at the keV level over this range. 

We first study the cutoff dependence of the LCS gaps in 
Fig.|U For simplicity, we here include only the contri- 
bution. We find that the cutoff variation is small overall 
and for 2.0 fm" 1 < A < 2.5 fm -1 the gaps are approx- 
imately cutoff independent. Next, we check our results 
against the gaps obtained in Ref. [l6j] using the same 
pairing interaction: the part of the smooth Viowfc 
evolved from the Argonne Dig potential to A = 2.5 fm -1 
with n oxp = 6 (and the same Skyrme functional SLy4 in 
the particle- hole channel). The very good agreement is 
shown in Fig. O The gaps are within 0.5% for almost 
all tin isotopes, with a couple of them (N = 54 and 56) 
showing a larger (but still small) difference of 3%. This 
demonstrates that the truncation after the first HFB iter- 
ation introduces a negligible error for pairing gaps. In the 
following, all results are based on the RG-evolved N 3 LO 
potential of Ref. (25| with A = 2.0 fm -1 and n cxp = 4. 



FIG. 6: State-dependent canonical pairing gap A aa in 120 Sn 
versus canonical single-particle energy e a - In addition, we 
show the occupation factor (labeled V 2 ), the 

E 9 UIVS factor (lab eled UV ), and the (£ q tW)V3?»+T 
factor (labeled UVy/2j + 1) that enters the gap equation, 
Eq. ((2}. The y-axis for the different factors is dimension- 
less and the factors corresponding to the three partial-wave 
cases are the same on this scale. 
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In Fig. [5J we plot the state dependence of the canonical 
pairing gap A Qa versus the canonical single-particle en- ^a^a)V^ja + 1 factor shown in Fig. [6) the effect of 
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FIG. 7: LCS and average gaps in the lead isotopes with in- 
creasing partial-wave contributions to the pairing interaction. 
Experimental gaps based on the three-point mass formula are 
shown for comparison. 



ergy e a . The occupation factor (V 2 ) and the UV factors 
are included to show the Fermi level. As expected from 
the pairing matrix elements, we observe that the impact 
of the 3 Pi partial wave is attractive for some states (pair- 
ing gap is increased) and repulsive for others (pairing 
gap is decreased). This is due to the interplay between 
upper-lower spin-orbit configurations for spin-triplet con- 
tributions to pairing, discussed in Sect. IIV Al Given the 
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FIG. 8: Same as Fig. [7] but for tin isotopes. 



FIG. 10: Same as Fig. [7] but for calcium isotopes. 
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FIG. 9: Same as Fig. [Jj but for nickel isotopes. 



FIG. 11: Same as Fig. [7j but for oxygen isotopes. 



the 3 Pi partial wave on A aa in the gap equation, Eq. 
is driven by the sign of the 3 Pi contribution to the pairing 
matrix elements (aa \ (I-P12) Mowfc | //) j=o,t=i,t m =-i, 
where / denotes the LCS state close to the Fermi level 
where (J2 q U£V£)^/2j a + 1 is peaked. Following the 
analysis of Sect. IIV Al for the LCS state (a = f) or if 
a and / are the same upper or lower spin-orbit config- 
urations, the 3 Pi partial wave will decrease the state- 
dependent pairing gap A aa , while for different spin-orbit 
states a, the 3 Pi contribution will increase A aa . More- 
over, as expected from Figs. Q] and [21 the impact of higher 
partial waves beyond the 3 Pi channel is small in Fig. [S] 
and for T = 1 pairing gaps in other nuclei. 

We present a more global study of neutron-neutron 
pairing gaps by calculating the LCS and average gaps for 
even-even nuclei in the lead, tin, nickel, calcium and oxy- 
gen isotopes in Figs.lTlfTTl The theoretical gaps are shown 
with increasing partial-wave contributions to the pairing 
interaction and in comparison to the experimental gaps 



based on the three-point mass formula, Eq. ((6]), centered 
on the neighboring odd nuclei. Higher partial waves were 
also included in Ref. [26(. The additional contributions 
beyond the standard 1 So channel are dominated by the 
3 Pi partial wave. They lead to a decrease of the LCS 
gaps of up to 15% and can change the isotopic depen- 
dence at this level. On the other hand, because of the 
interplay of spin-orbit configurations for the spin-triplet 
3 P! channel and the resulting alternation of repulsive and 
attractive contributions to the state-dependent pairing 
gap in Fig. we find that the average gap changes little 
when additional partial waves beyond the channel 
are included. This shows that the impact of the 3 Pi con- 
tribution depends on the definition of theoretical gap. 
We therefore conclude that future studies should com- 
pare directly calculated odd-even mass differences to the 
experimental gaps. This has recently been carried out in 
Ref. using only the 1 So contribution to the pairing 
interaction. 
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FIG. 12: T = 0, J = 1 pairing matrix elements in the sd shell 
based on the same RG-evolved N 3 LO potential as used to 
calculate the T = 1, J = pairing matrix elements in Fig. [T] 
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FIG. 13: Same as Fig. [12] but for the pf shell, where partial 
waves up to the 3 G3 channel can contribute to T = 0, J = 1 
pairing matrix elements. 



VI. T = 0, J = 1 PAIRING MATRIX ELEMENTS 

We turn to the case of neutron-proton T — 0, J = 1 
pairing, with a focus on the relative pairing strengths 
compared to the T — 1, J = channel. In Figs. [T^l 
and [T31 we show the jj-coupled T = 0, J = 1 pairing 
matrix elements in the sd and pf shell. The matrix el- 
ements are in the HO basis with hco = lOMeV, but the 
general effects of partial waves are qualitatively similar 
for other values of tko or in the HF basis. While all 
T = 1, J = pairing matrix elements are attractive at 
the S-wave level, for T = pairing the S-wave contribu- 
tion is given by the spin-triplet 3 Si — 3 Di channel, which 
is influenced by the interplay between spin-orbit config- 



FIG. 14: Comparison of T = 1, J = versus T = 0, J = 1 
pairing matrix elements in the pf shell. Interactions that 
favor T = 1 pairing lie below the diagonal. The only point 
above the diagonal corresponds to the (pi/2) 2 — (P1/2) 2 matrix 
element. The pairing matrix elements and the symbols used 
are the same as in Figs. [2land [l3l 



urations discussed in Sect. IIV Al As a result, although 
the spin-triplet S-wave interaction is attractive, pairing 
matrix elements connecting different spin-orbit configu- 
rations, {j a >) 2 -{jb<) 2 and (j a <) 2 -(jb>) 2 , are repulsive. 
This can be seen clearly in the (d, d) , (p, p) , (/, /) and 
(p, f) groups at the 3 Si - 3 Di level in Figs. QJ and [H 

The next partial-wave contribution is due to the re- 
pulsive *Pi spin-singlet channel (in contrast to the spin- 
triplet 3 Pi for T = 1 pairing), so that it decreases the 
attractive pairing strength of all T = 0, J = 1 pairing 
matrix elements. Finally, Figs. [12] and [13] show smaller 
contributions due to the spin-triplet D waves, again with 
alternating attractive and repulsive character. 

We compare the pairing strengths in the pf shell by 
plotting the T = 1, J = versus T = 0, J = 1 pairing 
matrix elements in Fig. [TD This clearly demonstrates 
that nuclear forces in the pf shell favor T = 1 over T = 
pairing, except for the (pi/2) 2 — (pi/2) 2 matrix elements. 
This is in contrast with the expectation based on the rel- 
ative S-wave interactions in free space, which favor the 
T = over T = 1 channel, with a bound deuteron com- 
pared to the nearly-bound state in the channel. 



A. Spin-orbit suppression 

In addition to the effects of the repulsive : Pi chan- 
nel and of D waves, the suppression of T = pairing 
at the S-wave level in Fig. [H] is due to spin-orbit split- 
ting between single-particle states. The impact of spin- 
orbit splitting onT = pairing properties has also been 
pointed out in shell-model calculations using empirical 
interactions [H, [36|. In our analysis, the significance of 
spin-orbit splitting lies in the choice of the j ' j-coupled ba- 
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FIG. 15: Ratio of the diagonal (luj,) 3 - (i„j a ) 2 pairing matrix 
elements in the T = 0, J = 1 over the T = 1, J = channel. 
Results are shown including only S waves as a function of the 
total angular momentum j up to the pf shell. The dotted 
line connects the p/-shell matrix elements. The solid line is 
the ratio R of Eq. for S-wave contact interactions with 
identical strengths. 



FIG. 16: Absolute value of the ratio of T = 0, J = 1 over the 
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T = 1, J = pairing matrix elements (all (l a j a 
in the pf shell. The anomalously large ratio for the (P1/2Y - 
(P1/2) 2 matrix element beyond S waves is due to the smal 
T = 1, J = pairing matrix element in this case (see Fig. [2} 



sis and in our assumption that pairing matrix elements 
based on low-momentum interactions provide a first ap- 
proximation for pairing in a j-sheil [HI, [l6| • 

We discuss this first at the level of the 1 So and 3 Si 
partial waves. The corresponding pha se shifts [30l ] (and 
also low-momentum interactions |{J [l0( ) have very similar 
momentum dependences, because at low energies nuclear 
forces mainly differ by providing slightly more attraction 
in the T = channel to lead to a loosely bound deuteron. 
The similarity of the 1 So and 3 Si channels is also re- 
flected in large scattering lengths, ai Sg — — 23.768fm 
and ct3 Sl = 5.420 fm. At low energies, nuclear interac- 
tions are therefore close to Wigner's SU(4) limit, where 
S-wave contact interactions have identical strengths [37| • 

In order to study the ratio R of the T = 0, J = 1 over 
the T = 1, J = pairing strengths, we consider two S- 
wave contact interactions with the same T = 1 and T = 
strengths. The resulting pairing matrix elements can be 
calculated analytically (see, for example, Ref. [38|]) and 
the ratio R is given by a geometrical factor: 



R = 



(njj njj I (1 - P12) V3 Sl I n b lj n b lj)j =1:T =a 
(njj njj I (1 - Pi 2) Vi So I n b lj n b lj} J=0tT =i ' 



j 3 1 
1/2 -1/2 



j j 1 
1/2 1/2 -1 



j 3 
1/2 -1/2 



4j 2 + 4j + 3 
8jC? + 1) 



(7) 



where (...) are 3j symbols [28j]. We find R = 1 for j = 
1/2, P = 1/2 in the limit of large j, and already R ^ 3/5 



for j ^ 3/2 [39(. This geometrical factor is due to the 
spin-orbit suppression and provides a simple explanation 
why S-wave nuclear forces favor T = 1, J = over T 
li. ./ : I pairs, except in j = 1/2 orbitals. 

In Fig. [151 we compare the ratio R to the diagonal 
((„jj 2 — (l a j a ) 2 pairing matrix elements up to the pf 
shell. We find that the trend of Eq. ([7]) agrees nicely. As 
expected, the ratio is somewhat larger for low-momentum 
interactions, because the 3 Si part is more attractive than 
the 1 S partial wave. For all j > 3/2, the T = 1, J = 
channel is favored. Finally, we include all partial-wave 
contributions and show in Fig. [TC] the absolute value of 
the ratio of T = 0, J = 1 over the T = 1, J = pairing 
matrix elements in the pf shell. This demonstrates that 
nuclear forces favor T = 1, J = over T = 0, J = 1 
pairing, except in low-j orbitals. 

The suppression R — 1 /2 in the limit of large j can be 
interpreted as a distribution of the 3 Si strength on more 
configurations, compared to the case, because there 
are only two states for J = 0: (j>) 2 and (j<) 2 ; compared 
to four states for J = 1: (j>) 2 , (j<) 2 , (j>i<) and (j < j>). 
In addition, a semiclassical picture for the suppression of 
T = pairing has been developed in Ref. [401 ]. 



VII. CONCLUSIONS AND OUTLOOK 

We have studied in detail partial-wave contributions 
of nuclear forces to pairing in nuclei, at the level of the 
pairing matrix elements and for neutron-neutron pair- 
ing gaps in semi-magic isotopic chains. For T = 1, J = 
pairing, the repulsive 3 Pi channel decreases the LCS gaps 
by up to 15% compared to the gaps obtained from the 
standard contribution, while the changes in the av- 
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erage gap were found to be small. The latter is due to 
the alternation of repulsive and attractive contributions 
for spin-triplet nuclear forces between paired nucleons. 
While we have focused on neutron-neutron pairing gaps, 
our conclusions equally apply to proton-proton gaps. 

We expect the 3 Pi effects will be more important for 
phenomena involving state-dependent pairing gaps, such 
as particle-transfer reactions. From the differences be- 
tween the LCS and average gaps beyond the 1 So channel, 
we conclude that future studies should compare directly 
calculated odd-even mass differences to the experimen- 
tal gaps (see also Ref. Finally, important future 
work for pairing gaps is to include the effects of three- 
nucleon forces (for the impact in neutron matter, see 
Ref. [14in . many-body contributions or induced interac- 
tions [19l 12011 . and center-of-mass corrections to pairing 
interactions [261 ]. 

Based on the comparison of pairing matrix elements, 
we have shown that nuclear forces favor T = 1 , J = over 
T = 0, J = 1 pairing, except in low-j orbitals. This is 
in contrast to the relative S-wave strengths in free space. 
We have traced the suppression of T = pairing to two 
origins. First is the spin-orbit splitting, which results 
in the 3 Si strength being distributed on more configura- 
tions, compared to the 1 So strength in the case of T = 1 
pairing. In the SU(4) limit where both relative S-wave 
interactions are equally attractive, this leads to a relative 



suppression of the pairing interaction for large- j orbitals 
by a factor 1/2, in favor of the T = 1, J = case. The 
deuteron 3 Si channel would have to be twice as attractive 
as to overcome this spin-orbit suppression. Second, 
T = pairing is weakened by the additional repulsive : Pi 
channel and by the effects of D waves, more strongly com- 
pared to the effects of higher partial waves in the T = 1 
pairing case. Our findings render the standard free-space 
motivation for T — pairing at best incomplete. 
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